GENERALIZED YANG-BAXTER EQUATIONS AND BRAIDING 

QUANTUM GATES 
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Abstract. Solutions to the Yang-Baxter equation - an important equation 
in mathematics and physics - and their afforded braid group representations 
have applications in fields such as knot theory, statistical mechanics, and, most 
recently, quantum information science. In particular, unitary representations 
of the braid group are desired because they generate braiding quantum gates. 
These are actively studied in the ongoing research into topological quantum 
computing. A generalized Yang-Baxter equation was proposed a few years ago 
by Eric Rowell et al. By finding solutions to the generalized Yang-Baxter equa- 
tion, we obtain new unitary braid group representations. Our representations 
give rise to braiding quantum gates and thus have the potential to aid in the 
construction of useful quantum computers. 



1. Introduction 

In 1944, L. Onsager published his now-famous solution of the Ising model in 
statistical mechanics. He employed in his solution a clever relation called a star- 
triangle transformation |PAj . better known today as the Yang-Baxter equation. 
Although Onsager never fully realized the significance of the star-triangle trans- 
formation, R. J. Baxter did, and he demonstrated its importance in the exact so- 
lutions of other statistical mechanical models. C. N. Yang, working independently 
to solve one- dimensional many-body problems with delta-function interactions, 
came across the same transformation. The term Yang-Baxter equation was coined 
by L. Faddeev in the late 1970s in honor of these two physicists. Now the Yang- 
Baxter equation is an important equation in mathematics and physics, connected 
to a variety of fields such as statistical mechanics, quantum field theory, quantum 
topology, quantum groups, and, most recently, quantum information science. 

The Yang-Baxter equation in dimension d (some would say in dimension d 2 ) is a 
matrix equation for an invertible complex matrix R = (R^),i,j, k, I = 1, 2, • • • , d. 
The easiest way to write down the Yang-Baxter equation is to use the language of 
linear operators between vector spaces. Let V be a <i-dimensional complex vector 
space with a chosen basis {e^}, % = 1, 2, • • • , d. The matrix R defines an invertible 
operator R : V ®V — > V ®V by R(ei®ej) = Ylki=i^ij e k® e hhj = 1> --,d, where 
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{e a <8> e^}, a, b = 1, .., d, is a basis of V <8> V. By abusing notation, we have denoted 
the operator associated to i? also by R. Letting I v be the identity operator on V, 
we form two operators R® Iy and Iy ® R from V <g> V <8> U to itself. Then the 
Yang-Baxter equation, written as an equation for linear operators, is the following: 

(YBE) (R ® Jy)(/y ® R)(R ® ly) = (ly ®R)(R® Iy)(/y ® i?). 

When we use the basis {e a ® e& <8> e c }, a, 6, c = 1, d, the above equation for 
linear operators becomes an equation for the matrix R. This matrix equation 
for R is usually called the Yang-Baxter equation, and a solution is called an R- 
matrix. The two equations are in one-to-one correspondence when bases of the 
involved vector spaces are fixed, so we will speak of both as the Yang-Baxter 
equation with the understanding that some basis has been chosen for the operator 
equation. Strictly speaking, the Yang-Baxter equation here is the braided version 
of the constant quantum Yang-Baxter equation, but we refer to it simply as the 
Yang-Baxter equation for brevity. 

The Yang-Baxter equation consists of a set of polynomial equations for the 
entries {Rfj} of the matrix R. Explicitly, for any choice of two ordered triples 
(x, y, z) and (u, v, w) from {1, 2, • • • , d}, we have 

f~Y"R'p 1 '\ \r~*d P a b f? cz f? x H ST^-d r>np J3xm T>yz 

\ / Z_/a,6,e=l ^uv^bw^ac Z—/m,n,p=l ^vw^un - rt mp - 

One application of unitary i?-matrices is to quantum information science. A 
unitary i?-matrix leads to a unitary representation of the braid group, and the 
resulting unitary matrices associated to braids can be used to process quantum 
information [NSSFDJ. Inspired by such an application, E. Rowell et al proposed 
a generalized version of the Yang-Baxter equation in [RZWGJ. Unitary solutions 
to the generalized Yang-Baxter equation sometimes also afford braid group repre- 
sentations and therefore can be used for quantum information processing. An 
8x8 solution to a generalized Yang-Baxter equation is used to generate the 
Greenberger-Horne-Zeillinger states [RZWG] . 

Solutions to the Yang-Baxter equation are difficult to find. This can be seen 
by counting the number of variables and the number of equations. If the vector 
space V is of dimension d, then R has d 4 entries, so there are d 4 unknowns. The 
Yang-Baxter equation consists of d 6 cubic polynomial equations for the d 4 variables 
{Rfj}- For d — 1, the equation is automatically satisfied by any nonzero complex 
number. For d = 2, the Yang-Baxter equation consists of 64 cubic homogeneous 
polynomial equations for 16 complex unknowns. This is the only case for which 
the Yang-Baxter equation has been solved completely, albeit with the help of 
a computer. The general solutions are found in [Hj, and based on this work, 
the unitary ones are classified in [Dj. The case of d = 3 consists of 729 cubic 
polynomial equations for 81 unknowns. No wonder it has never been completely 
solved! 



3 



Furthermore, new solutions to generalized Yang-Baxter equations are even more 
challenging to discover. So far, only two essentially new solutions have been found: 
the 8x8 unitary solution given in [RZWGJ and a unitary solution found by 
Prof. Rowell based on an unpublished work of Prof. Goldschmidt and Prof. Jones 
(cf. |GHP4 lR]). We will refer to the latter as the Rowell solution. In this paper, we 
generalize the Rowell solution to three families of unitary solutions, each of which 
is parameterized by the upper semicircle in the complex plane. Our solutions lead 
to new unitary braid group representations that generate braiding quantum gates. 

Our solutions were found with the help of the computer package Maple. Mod- 
ifying a short program written by Prof. Rowell, we found the solutions up to an 
error of order 10 -11 . Then we proved algebraically that they are indeed solutions. 

The contents of the paper are as follows. In Section [2], we review the generalized 
Yang-Baxter equations. Section [3] contains our main result: certain solutions to a 
generalized Yang-Baxter equation. In Section HI we show that all of our solutions 
in Section [3] give rise to unitary representations of the braid group. In Section O 
we list some open problems and directions for future research. 

2. Generalized Yang-Baxter Equations 

The Yang-Baxter equation (YBE) is indexed by a single natural number d, the 
dimension of the vector space V. The generalized Yang-Baxter equation (gYBE) 
proposed in [RZWGj is indexed by two more natural numbers, m and I. For 
convenience, we will continue to use operator language while keeping in mind 
that, by the choice of a basis of the vector space, the gYBE is a matrix equation. 

Definition 2.1. Let V be a complex vector space of dimension d. The (d,m,l)- 
gYBE is an equation for an invertible operator R : V® m — > V® m such that 

(gYBE) (R ® lf)(lf ®R){R® If) = {if ®R){R® lf){lf <8> R), 

where d, m, and I are natural numbers and if is the identity operator on V® 1 . 
Any matrix solution to the (d, m, Z)-gYBE is called a (d, m, /)-i?-matrix. 

Note that if m = 2 and 1 = 1, the gYBE reduces to the usual YBE. Generally, 
however, the gYBE is harder to solve than the YBE, with the exception of a couple 
of basic cases. When d = 1, R is just a scalar, so any non-zero complex number is 
a solution. When m = I = 1, the gYBE becomes the equation R 2 <g> R = R® R 2 , 
where R is an invertible operator on V. R 2 ®R = R®R 2 is the same as R®I = I®R 
because R is invertible. It follows that R = \Iy for some nonzero scalar A, if 
m = I = 1. 

For the application to quantum information science, we will focus on d = 2 
because, when d = 2, V is isomorphic to C 2 , the so-called qubit state space. The 
(2, 2, l)-gYBE is the YBE in dimension 2, so the first non-trivial gYBE for qubits 
is the (2, 3, l)-gYBE. 
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Definition 2.2. The (2, 3, l)-gYBE for R : (C 2 ) 03 (C 2 ) 03 is 

R1R2R1 = R2R1R21 

where Ri = R ® I 2 and R 2 = h®R act on {C 2 ) m . 

There are several different conventions regarding the tensor product of matrices. 
We will use the so-called Kronecker product: for two matrices A = (aij) mxn and 
B = {bki) p xq, A £g> B is the (rap x ng)-matrix obtained by replacing each entry 
of A by the block a^B. We will use /„ to denote the n x n identity matrix. When 
no confusion would result, we will simply write I n as I. 

(X 
Y 

Proposition 2.3. Suppose X is a Ax A matrix. Then X © X is a solution to the 
(2, 3, l)-gYBE if and only if X is a solution to the YBE in dimension 2. 

A proof is contained in the proof of Lemma 13.51 

As mentioned in the introduction, there currently exist only two essentially 
new solutions for any gYBE other than the (d, 1, l)-gYBE. Up to simple changes 
described in Prop. 12.41 below, the first solution is a (2, 3, 2)-i?-matrix in Section 
4.2 of [RZWGj : 
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We will refer to this solution as the X-shape solution as the non-zero entries 
form the shape of an X. 

The second solution is the Rowell solution, a (2, 3, l)-_R-matrix. Let ( = e 2n ^ 8 . 
Then, 



Rr 



Any solution to the gYBE, just as for the YBE, leads to many more solutions 
by the following proposition. 

Proposition 2.4. If R is a solution to the (d,m,l)-gYBE, then 

(1) Ai? is also a solution for any nonzero scalar A. 

(2) is also a solution. 
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(3) (Q- 



R ■ Q® m is also a solution, where Q is an invertible dx d matrix. 



The proof for the YBE in [D] works for the gYBE without any change. The 
proof of the third case is especially interesting. 

One of the reasons for interest in the YBE is that any solution leads to a matrix 
representation of the braid group. This is not generally true for solutions to the 
gYBE, a point that we discuss in Section |U though the X-shape solution R x and 
the Rowell solution do lead to braid group representations. 

3. Solutions to the (2,3,1)-Generalized Yang-Baxter Equation 

In this section, we focus on the gYBE for (d, m, /) = (2, 3, 1) and generalize the 
only known non-trivial (2, 3, l)-i?-matrix: the Rowell solution R^. Our strategy 
is to search for unitary (2, 3, l)-i?-matrices with the non-zero entries in the same 
positions as in R^. We found that the Rowell solution lies in a one-parameter 
family and discovered two more families of new unitary solutions. In Section HJ 
we show that these solutions give rise to unitary braid group representations. 

The Rowell solution can be rewritten as 
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By Prop. 



the Rowell solution multiplied by ( is also a solution. 



Definition 3.1. Two (2,3, l)-i?-matrices R and R' are 

(1) locally conjugate if there is an invertible 2x2 matrix Q such that R = 
(Q- 1 ®Q- 1 ®Q- 1 )R'(Q®Q®Q). 

(2) equivalent if R and R' are related by a sequence of applications of the 
three relations in Prop. 12.41 (i.e., R can obtained from R' by multiplying 
a non-zero scalar, by taking the inverse, by a local conjugation, or by an 
arbitrary combination of these three operations). 



1) (2 x 2)-diagonal if M 



such that Mij,i,j = 1,2, are all 



Definition 3.2. A 4 x 4 matrix M is 

Xi M 12 
M 21 M 22/ 
diagonal 2x2 matrices. 
(2) (2 x 2)-diagonally unitary if it is (2 x 2)-diagonal and unitary, and M = 

•73 (^M ikf 2 ) SUC ^ eaC ^ ^ oc k i = 1,2, 3, 4, is unitary. 

Our main result is the following: 

Theorem 3.3. If an 8 x 8 unitary matrix solution R to the (2,3, l)-gYBE is of 
the form R = X © Y , where the 4x4 matrix X is (2 x 2) -diagonally unitary, then 
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'1) R is equivalent to an R(9) in one of the following three families for some 
< 9 < tt: 
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(2) Any two different (2,3,l)-R-matrices in the three families above are not 
equivalent to each other. 

(3) For each (2, 3, l)-R-matrix above, X is different from Y except when 9 = n 
in the third family. Therefore, neither X nor Y is a solution to the YBE 
unless for 9 = tt in the third family. 

For < 9 < tt, e %e is the upper semicircle in the complex plane. Therefore, 
solutions in each family correspond to points in the upper semicircle. For example, 
the Rowell solution corresponds to the point i in the complex plane; i.e., the Rowell 
solution is equivalent to R(tt/2) in the first family. When 9 = tt in the third family, 
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is a YBE solution that is locally conjugate to the i?-matrix corresponding to the 
Bell states [D] [KL]. 



3.1. Proof of Main Theorem 13.31 We start the proof of Thm. lSTBl by presenting 
several lemmas. 

In this subsection, we use R,X,Y, A, B,C, D,Yi,i = 1,2,3,4, and a,f3,cj,j,8 
to denote the following matrices and matrix entries. The 8x8 matrix R is X © 
Y for some 4x4 matrices X and Y. The matrices X and Y are written as 



7 0\ 

q ^ J . We also denote the identity matrix 7 2 simply by 7. 



Lemma 3.4. If X is (2x2) -diagonally unitary, then C = —DB^A 



-ja 
-Sup 



Proof. XX* = I A implies that CA* + 7)77+ = 0. Thus, C = -DB\A^)-\ If A is 
unitary, then (A^)^ 1 = A. A simple computation yields C in terms of a, j3, u>, 7, <5. 

□ 

Lemma 3.5. in 8 x 8 matrix R = X ®Y is a solution to the (2, 3, Y)-gYBE if 
and only if the following equations are satisfied: 

(1) (A ® 7)X(A <g> I) + (B ® 7)Y(C <g> /) = >/2 • X(A ® 7)X 

(2) (A ® 7)X(77 <g> 7) + (77 ® 7)Y(L> <g> /) = y/2 ■ X(B ® 7)Y 

(3) (C <g> 7)X(A <g> J) + (D <g> 7)Y(C <g> 7) = y/2 ■ Y{C ® 7)X 

(4) (C <g> 7)X(77 <g> 7) + (D ® 7)Y(L> <g> 7) = • Y (D ® 7)y 

(5) (y <g> 7)x(y ® 7) + (y 2 <g> 7)y(y 3 ® 7) = • x(y ® 7)x 

(6) (y ® 7)x(y 2 ® 7) + (y, ® 7)y(y 4 ® 7) = V2 ■ x(y 2 ® i)y 

(7) (y 3 ® 7)x(y ® 7) + (y 4 ® 7)y(y 3 ® i) = V2 • y(y 3 ® 7)x 

(8) (y 3 <g> 7)x(y 2 ® 7) + (y 4 ® 7)y(y 4 ® 7) = V2 ■ y(y 4 ® i)y 

Proof. Recall that the (2, 3, l)-gYBE is RiR 2 Ri = R2R1R2 for Ri = R ® 7 = 
(f ^ Y ® 7^ an< ^ Fi2 = I ® R = ■ The left-hand side is 

X®I \(R 0\(X®I \ _ f(X ® 7)77(X ® 7) 

y®7 y )\ v o77 y )\ v o y®7 y )~\ v (y ® 7)77(y ® i) / 

while the right-hand side is 

R 0\(X®I \fR 0\_(R(X®I)R \ 

RJ \ Y®l)\0 RJ ~ \ R(Y®I)RJ- 

Hence the (2, 3, l)-gYBE for 77 = X © Y is equivalent to 

(1) (X ® 7)77(X ® 7) = 77(X ® 7)77, 

(2) (Y ® 7)77(y ® 7) = R(Y ® 7)77. 

Substituting X ® I = ^ ® ^ 7) ® 7^) anc ^ ^ = y^ * n ^° ^ e ^ rs ^ ec L ua ~ 
tion, we get 



(A®I B®l\(X $\(A®I B®I\_ rrfX 0\(A®I B®I\[X 0\ 
\C®I D®l)\0 YJ\C®I D®l) V l v Y)\C®I D®l)\0 YJ 
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Comparing the blocks, we obtain the first 4 equations. By replacing A, B, C, D 
with Yi, i — 1,2, 3, 4, respectively, we obtain the second 4 equations. 

□ 

Lemma 3.6. Suppose an 8x8 matrix R = X@Y is a solution to the (2, 3, l)-gYBE 
such that X is (2 x 2) -diagonally unitary. Then, 

(1) Equation (1) of Lemma \3.5\ is equivalent to 

Mn = -7^-^-/)=(£ W , (1+ L_ w) )- 



(b) Y 2 = -a(38uB(A + cjD - ul) -- > _ ,., , , 



^(1-7-w) 
-a/3 2 
/3(l + w7-w) 

a/3 2 5 2 w 2 7 ' ' 



(c) F 3 = a/375t J DA(A + wD-w/)=^ 

( d) n = sm^ - wD + „») = - t) t _ J + J . 

(2) Equation (2) of Lemma \3.5\ is equivalent to 

(a) (8 - l)u 2 + (1 + 5 2 - 5y + 7 )w - 2 7 = 0, 

(b) 5(2 — — 7) = — a) — 7 + I + W7 — 7 2 + cjj, 

(c) (5-l) 7 = 5 2 -l + a;(l-5), 

(d) 5(2cj + ui — 7) = w — I + 7 + W7 + W7 — 7 2 . 

(3) X zs a solution to the YBE in dimension 2 if and only if X = Y . 

(4) Y is unitary if and only if 

(a) w + a; + 7 + 7 = + £7 + 2, 

(b) w + w + 5 + 5 = u5 + u5 + 2, 

(c) 1 + lo + lo + wy = 7 + 7 + + ^7, 

(d) 2 + w5 = 5 + 8 + u5, 

(e) u; + u + 7 + 7 + ^7 + ^7 = 4 + cu 2 + u 2 , 

(f ) 5 + 5 + w5 + w5 = 2 + w + u. 

(5) y zs unitary if and only if to, 7, 5 fall into one of the following three cate- 
gories: 

(a) to = 7 = ±i, 8—1, 

(b) w = 5 = ±i, 7 = 1, 

(c) to = 7 = 8 = 1. 



Proof. For (1), by equation (1) of Lemma | 

y = (£t <g> /) [\/2 ■ X(A g> I) X — (A<S> I)X(A ® /)] (C f <g> /) . 

Substituting A®J=g ^^^(f ^).^=(T J -/&)/)■ 
and X = -j= into the equation, we obtain 



j_ / - 7w5C + 7 A -a8/3uAB - a5/3BD + ad/3B\ 

V2 \-^ a (3CA - jaficuDC + ^/af3cuC -5cuCB - 5D 2 + 5uD J ' 

The formulas follow from replacing C with — DB^A and simplifying. 

For (2), we use equation (2) of Lemma [331 

For (3), X satisfies the YBE in dimension 2 if and only if X satisfies the first 4 
equations in Lemma [3.61 when Y is replaced by X. Now (3) follows. 
For (4), YY* = I4 is the same as 

(1) Y 1 Y}+Y 2 Y} = 2I, 

(2) y x y£ + Y 2 Yj = 0, 

(3) Y 3 Y} + Yjr} = 0, 

(4) Y 3 Y* + Y 4 Y} = 21. 

Substituting Yj, % — 1, 2, 3, 4, from above gives the equations. 

For (5), by (d) above, we have 2 — 5 — 5 = u{5 — 5). It follows that oj{6 — 5) is a 
real number. Let u = x + iy,5 = a + ib. Since (2 + iy)(—2ib) is real, its imaginary 
part is —2x6 = 0. If a; = 0, then 00 = ±1 If 6 = 0, then 2 — 5 — 5 = 0; hence 5 = 1. 

Case 1: u — i 

The equations for Y to be unitary become 

(1) (l-i) 7 +(l + i)7 = 2, 

(2) (l-i)5 + (l+i)5 = 2. 

Thus, the real part of (1 — 2)7 is 1. When we let 7 = c + cfi, c — d = 1. Since 
c 2 + d 2 = 1, it follows that either c = 1, d — or c = 0, d = 1; i.e., 7 = i or 7 = 1. 
Similarly, 5 — i or 5 = 1. 

If 5 = z, then equation (c) of (2) implies that 7 = 5 + 1 — i — 1. Therefore, 
u — 5 — i, 7 = 1. 

If 5 = 1, substituting into (a) of (2), we get u = 7 = i. Therefore, u = 7 = 
z,5 = l. 

Case 2: u = — i 

This case is the complex conjugate of Case 1. Therefore, uj = 5 = — i, 7 = 1 or 
u; = 7 = — i, 5 = 1. 
Case 3: 5 = 1 

By (a) of (2), we have u = 7. Then the equations for Y to be unitary become 
2(u; + a)) = u 2 + Co 2 + 2. Hence u + w = or 2. It follows that co — ±z or u — 1. 
Therefore, cj = 7 = ±i, 5 = loru; = 7 = 5 = l. 

□ 

Lemma 3.7. Le£ i? = X © Y, u>z£/i X and Y being the conjugates of X and Y 
as follows: X = g °) X g and Y = g _ ^ Y g J^J . Lei 

A, £?, C, D, Y, 2 = 1,2, 3, 4, denote the corresponding matrices in R. 

If R is a solution to the (2,3,l)-gYBE and X is (2 x 2) -diagonally unitary, 
then R is also a solution to the (2,3,l)-gYBE with B being the identity matrix 
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I. Conversely, if R is a solution to the (2,3,l)-gYBE such that B = I and X 
is (2 x 2) -diagonally unitary, and B = is an arbitrary unitary diagonal 

matrix, then R = X@Y , which is similarly obtain from R and B, is also a solution. 
Proof. Direct computation shows that 



X 



J_ ( A 1 
^2 {-DA D 



Y 



-( Yl 

V2 \a/3BY 3 



a/3B*Y 2 
Y 4 



Therefore, 

A = A, B = I, C = -DA, D = D, 
Y x = Y 2 , Y 2 = aPB^Y 2 , Y 3 = a(3BY 3 , Y 4 = Y 4 . 

Replacing all matrices in Lemma 13.51 by the R matrices using the above identi- 
ties, we find that R is a solution to the (2, 3, l)-gYBE if and only if 



(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 



A 

A< 



y/2-X{A®I)X 



I)X(A® I) + Y{{-DA) ® I) 
I)X + Y(D <g) /) = y/2 ■ XY 
{(-DA) <g> I)X(A ®I) + (D® I)Y((-DA) <g> /) = y/2 ■ Y((-DA) g 
{(-DA) <g> I)X + (D® I)Y(D <g> I) = V2 • Y(D g> I)Y 
[Yi ® I)X(Y 1 <g> /) + ((afiB^Y 2 ) <8> I)Y((af3BY 3 ) ® I) = V2 ■ <g 
(Y 1 ®I)X((a^B^Y 2 )®I) + ((a(5B^Y 2 )®I)Y(Y A ®I) = V2-X((a^Y 2 ) 
(a(5BY 3 )®I)X(Y l ®I) + (Y 4 ®I)Y((a(5BY 3 )®I) = V2-Y((af3BY 3 ) 
(a(3BY 3 ) ® I)X((a(5B^Y 2 ) ® I) + (Y 4 ® I)Y(Y 4 ® /) = y/2 ■ Y(Y 4 (g 

It suffices to show that the above 8 equations are the same as the 8 equations of 
Lemma 13.51 

The 8 cases are similar, so we show only the first one in detail. Substituting X 
and Y into the first equation, we have 



)I)X 

>I)X 
®I)Y 
%I)X 
I)Y 



(A® I) 



I 
B 



X 



I 

B^ 



*-\[ b) x Co B t 



(A ® 7) 







I 

a(5B 



Y 



I 

a@B* 



((-DA)® I) 



Multiplying the equation from the left by ^ and from the right by 





1 






G 


Bt J' 


5 


5 



another, we obtain 
(A®I)X(A®I) 



(0 ^/) r (o Q g J )(-^( 8 j) = ^(A®Dx 
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Since 



a(B g) J) and 







(al 




o 


) 


= a (o 


° ) 

01 



I \ _ _ (al 

a (B^ (g) J), we have 

(A © J)X(A <g> J) + (£ <g> I)Y(-B ] DA g> J) = v / 2X(A © J)X, 
which is equation (1). 



□ 



In particular, Lemma 13.71 shows that any solution R = X 
diagonally unitary can be recovered from a solution with £> ; 
and 0. In the unitary case we have the following: 



■ Y with X (2 x 2)- 
J and the values a 



Lemma 3.8. The matrices R = X © Y for a 
7 = 1, uo — 5 — i; and a = f3 = u = ^ = 5 = 
solutions, respectively, to the (2,3,l)-gYBE: 

(1) 



R 



- — 6 — l,oo — j — i; a — — 
1 are the following three unitary 



(2) 



R 



(3) 





/l 





1 


o\ 




f » 





1 


o\ 


1 







i 





1 


i 





1 





-1 








— i 





i 







— z 





1 









V 





— i 





V 




v° 


i 





i 


/ 






/ 1 





1 






(i 





1 






i 







/ 





1 


1 





1 





-1 




= 71 




-1 





1 





© -P 
\/2 


i 





i 











V° 


1 









V° 


1 





l ) 






^ 1 





1 


o\ 




A 





-1 


o\ 


i 







1 





i 


1 





1 





-1 






-1 





1 





© -p 


i 





1 









V 





-1 





V 




v° 


1 





1 


/ 


eluded the 


solutions for a = 






5 = 


= 1, 


UO 



7 



and 



a = (3 = j = l,u 
(2), respectively. 



5 = i, which are merely the conjugates of solutions (1) and 



in 



Proof. It suffices to show that the three matrices satisfy the equations (1) — 
Lemma 13.51 

The verification is a long matrix algebra computation. This can be done using 
a software package such as Maple. To do the computation algebraically, we notice 
that the verification is equivalent to checking 8 sets of matrix equations. We have 
done both the Maple verification and the algebraic computation. □ 



Lemma 3.9. The following are true: 
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(1) For each (2, 3, l)-R-matrix R = X © Y in Lemma \3. 81 X and Y are con- 
jugate to each other. Their common eigenvalues are {e~^, e~i2 , e~, e^"} 
for the first family, {e - "?, — e~^, e"^, e^} /or the second family, and {e - ^, e" 
for the third family. It follows that any two different (2,3, l)-R-matrices 
above are not conjugate to each other. 



(2) For any a, (3 such that \a\ 



1, the matrix R(a } (3) in each family 



below is a unitary solution to the (2,3, l)-gYBE: 
(a) 







i 





a 




1 







i 





/3 


71 




-ia 





% 







V 





-if3 





V 







i 





/3 


\ 


i 







1 





-a/3 2 






-i/3 





1 







V 





iaf3 2 





i ) 



(b) 



(c) 







/ 1 





a 






(i 





/3 


\ 


1 







i 








1 





i 





-a/3 2 


V2 


—a 





1 





v/2 







i 









v° 


£ 









V° 


a/3 2 





1 / 




/ 1 





a 






A 





-/3 





1 







1 





/3 


1 





1 





-a/3 2 


72 




-a 





1 





© -= 
V2 







1 







V 





-p 











a/3 2 





1 



/ 

(3) In each family above, R(a,(3) is locally equivalent to R(a',(3') if and only 

<i rv rv' 



Proof. For (1 
where a x = 



Y 

1 

1 



P^XP by P equal to 




ia, 

1 











and 



/ 
/ 



and cr 



-1 



are the Pauli matrices. 



The characteristic polynomials of X for the three cases are (A 2 -exA + *) 2 , 
((* - 73) 2 - |)((A - ^? + §), and ((A - ^) 2 + |) 2 . 
The solutions in (2) follow from Lemma 13.71 

a h 



For (3), suppose Q 



c d 



is an invertible 2x2 matrix such that 



(Q © Q © Q)R{a', /3') = i?(a', /3')(Q ©Q©Q). 

This equation is the same as 

(1) a{Q © Q)X(a', f3') = aX(a, (3){Q © Q), 

(2) 6(g © Q)Y(a', f3') = bX{a, (3)(Q © Q), 
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(3) c(Q ® Q)X(a', p') = cY(a, p)(Q ® Q), 

(4) d(Q ® Q)Y(a', 0') = dY{a, P)(Q ® Q). 

If neither a nor b is 0, then (1) and (2) imply X(a',(3') = Y(a',j3'). This contra- 
diction implies that one of a, b is 0. Since Q is invertible, therefore, exactly one of 
a, b is 0. Similarly, using (3) and (4), we know that exactly one of c, d is 0. Since Q 
is invertible, there are exactly two cases: b = c = and a 7^ 0, <i 7^ 0, ot a = d = 
and 6 ^ 0, c ^ 0. 

IfQ= (o °)' then 

(1) aQA(o^,p') = aA(a,P)Q, 

(2) aQB{a',P') = dB(a,p)Q, 

(3) dQC{a',P') = aC{a,(3)Q, 

(4) dQD{a',(3') = dD{a,(3)Q. 

Noting that Q, A, B,C, D all commute with one another, we have A (a', /?') = 
A(a,P),B(a/,p') = %B(a,P) and C(a',P') = %C(a,P),D(a!,P') = D(a,P). 
Therefore, R(a, (5) is locally equivalent to R(a', (3') by a diagonal Q if and only if 

I = §L 

a a' ' 

IfQ=(° J^then 

(1) 6gr 3 («',/3 , ) = c J B(a,/3)Q, 

(2) QYt(a',P') = A(a,/3)Q, 

(3) QFiK,/?') = D(a,P)Q, 

(4) C gy 2 (o ! , ,/3 / ) =bC(a,P)Q. 

For the first family of solutions, this is impossible. For the second and third 
families, the equations lead to the same conclusions as for the diagonal Q. 

□ 

Completion of the Proof of Theorem 13.31 If R is a (2, 3, l)-i?-matrix, by 
multiplying by a scalar, we may assume that the (1, l)-entry is 4j. Therefore, it is 
one of the R(a, P) in Lemma I'd. 91 By Lemma T3.9[ any two R-matrices in the same 
family are locally conjugate if and only if ^ are the same. Setting e %e = -, we know 
that all local conjugation classes up to a scalar are represented by < 9 < 2tt. 
Complex conjugation reduces the equivalence classes to < 9 < ir. 

□ 

4. Braiding Quantum Gates 

In quantum information science, information is encoded in quantum states and 
processed by unitary evolutions. In topological quantum computation, unitary 
evolutions are effected by braiding anyons, 2-dimensional quasi-particles |NSSFD] . 
Braiding anyons gives rise to unitary representations of the braid group. All of our 
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solutions to the (2, 3, l)-gYBE lead to unitary representations of the braid group 
and therefore can be used in quantum information science. 

4.1. Far Commutativity. Recall that the n-strand braid group B n has a pre- 
sentation with n — 1 generators <7i, (72, cr n -i, as well as two defining relations: 

(1) <7i<7 i+ i<7i = a i+1 (Ji(Ti + i for 1 < % < n - 2 

(2) <7j<7j = UjCTi for |i — j\ > 2. 

The first is commonly called the braid relation and the second far commutativity. 

A unitary i?-matrix naturally gives rise to a unitary representation of the n- 
strand braid group B n on V® n as follows. Set 

R . = 0j R(g) J® 

Then the assignment pR(o~i) = Ri is a unitary matrix representation of f? n . The 
braid relation 

Pi?(o"i)PR(o"i+l)PR(o"i) = PR(o"i+l)Pi?(o"i)p j R(o"i+l) 

follows exactly from the YBE for R. The far commutativity relation PR{o~i)pR{aj) = 
Pn((jj)pji{(Ti) j\> 1 holds because Pij(crj) and pn(crj) act nontrivially on dis- 

joint tensor factors of \/® n . 

If i? is a solution to a (d, to, Z)-gYBE, we can perform a similar assignment: for 
B n , set 

R a . = J®^" 1 ) ® R^I® 1 ^-^. 

Then the braid relation holds because of the gYBE for R. But the far commuta- 
tivity relation is not necessarily satisfied when > 1. However, if 21 > to, then 
it is still true that pR(o~i) and PR(o~j) act nontrivially on disjoint tensor factors of 
y®(m+(n-2)j) anc j {.j^ f ar comm utativity condition is satisfied. It follows that 
the Y-shape (2, 3, 2)-i?-matrix R x leads to unitary braid group representations. 
Otherwise, even when \i — j\ > 1, the two matrices Pi?(a"j) and PR^Uj) might not 
commute with each other. 

Definition 4.1. Suppose R is a (d, to, Z)-i?-matrix. Then i? satisfies the far com- 
mutativity condition if 

RaiRfTj RujRd\ 

in the braid group -Bq-— i)?+2 for every j such that j > 2 and (j — 1)/ < to. 

If a (d, to, /)-i?-matrix satisfies the far commutativity condition, then it yields 
a representation of each n-strand braid group B n on \/^ m +( n - 2 W. In particular, 
any (d, m, /)-i?-matrix with 21 > m leads to a representation of the braid group. 

Proposition 4.2. Let R = X®Y be a (2, 3, l)-R-matrix such that X = 

andY = ■+= ^ ^ for 2 x 2 matrices A, B, C, D, Y h i = 1,2, 3, 4. 
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(1) If A, B, C, D, Yi, i = 1,2, 3, 4, are all diagonal, then R ai R a3 = R a3 R ai . 

(2) // one of A, B, C, D, Y iy i = 1,2, 3, 4, is not diagonal, then R ai R a3 = R a3 R ai 
implies X = Y . 

Proof. The far commutativity condition is equivalent to 

(1) (X ® h)(I ®R) = (I®R)(X® h), 

(2) (V ® h){I ®R) = (I®R)(Y® I A ). 

These equations are the same as 

(1) {A®h)R = R(A®h), 

(2) {B®h)R = R{B®h), 

(3) {C®h)R = R(C®h), 

(4) (D®h)R = R(D®I A ), 

(5) (Yi ® I 4 )R = R(Yi ® J 4 ), % = 1, 2, 3, 4. 

Let A = ( ° 12 | . Then the first equation is 

\ a 21 ^22 y 

(1) anX = a\\X 

(2) a 12 F = a 12 X 

(3) a 2i X = a 2 iF 

(4) a 22 Y = a 22 Y. 

Now the proposition follows. 

□ 

Corollary 4.3. All (2,3,l)-R-matrices in Thm. \3.3\ satisfy the far commuta- 
tivity condition and thus yield unitary representations of the braid group B n on 

/(£2\(g)(n+l) 

4.2. i?-braiding Quantum Gates. 

Definition 4.4. Let R be a unitary (2, 3, l)-i?-matrix that leads to a representa- 
tion pr of the braid group. A unitary matrix Ui G U(2 n+1 ) is an i?-braiding gate 
if Ui — A • Pfi(cTj) for some braid generator o~i and non-zero scalar A. 

A unitary matrix U is an i?-braiding quantum circuit if there exists a braid h 
such that U = \ ■ pn(b) for some non-zero scalar A. 

Quantum circuits are compositions of quantum gates and are used to perform 
computations. It follows that each (2, 3, 1)-R matrix in Thm. 13.31 leads to (n + 1)- 
qubit braiding quantum circuits from its associated representation of the n-strand 
braid group B n for all n > 1. 

5. Open Questions 

Generalized Yang-Baxter equations provide a new method for finding represen- 
tations of the braid group. These representations have applications in a variety of 
fields; in particular, the resulting braiding quantum circuits are actively studied 
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in quantum information science [NSSFD] |KL] |RZ WG] . In this section, we list a 
few of the many open questions and directions for further research. 

(1) How can we find all unitary (2, 3, l)-i?-matrices R = X @Y such that X, Y 
are 2x2 diagonal and unitary? Note that their 2x2 blocks A, B, C, D, Yi, i = 
1, 2, 3, 4, are not necessarily unitary. 

(2) How do we find all unitary (2, 3, 1)- or (2, 3, 2)-i?-matrices with the zero 
entries in the same positions as in the solution Rx? 

(3) Which quantum circuits can be realized by the braiding quantum circuits 
resulting from our (2, 3, l)-i?-matrices? Mathematically, the question is to 
find the images of the afforded braid group representations. 

(4) Which entangled states can be generated from product states using our 
braiding quantum circuits? 

(5) For the three families of solutions in Thm 13.31 do any two (2, 3, l)-R- 
matrices R and R' from the same family lead to equivalent braid group 
representations? To determine whether the representations from R and Rl 
are equivalent, we need to find out whether there exists a single matrix P 
such that 

p-^R® I)P = R' ® I, P- 1 (I®R)P = I®R!. 

We hope to settle this question in the future. Note that the answer has no 
bearing on the application to braiding quantum gates. 

(6) Are our new representations of the braid group equivalent to known braid 
group representations obtained from other methods such as quantum groups? 

Ultimately we would like to know if there are real physical systems that would 
realize our braid group representations. Such physical systems would be quantum 
computers, powerful tools for exploring the quantum world and for bringing us 
new technologies of great benefit to society. 
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